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A disorder-dependent Gaussian variational approach is applied to the 
problem of a d dimensional polymer chain in a random medium (or poten- 
tial). Two classes of variational solutions are obtained. For d < 2, these 
two classes may be interpreted as domain and domain wall. The critical ex- 
ponent v describing the polymer width is v = j^z^j (domain solution) or 
v = (^4) (domain wall solution). The domain wall solution is equivalent to 
the (full) replica symmetry breaking variational result. For d > 2, we find 
v = No evidence of a phase transition is found for 2 < d < 4: one of the 
variational solutions suggests that the polymer chain breaks into Imry-Ma 
segments, whose probability distribution is calculated. For d > 4, the other 
variational solution undergoes a phase transition, which has some similarity 
with B. Derrida's random energy models. 
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I. INTRODUCTION 



Usually, in quenched disordered systems, extensive thermodynamical quantities, such as 
the free energy, are identified to their average over the disorder JIJ. This can be understood 
in the following way: one divides the macroscopic system into mesoscopic subsystems, each 
subsystem corresponding to a particular disorder configuration. For short-range forces, the 
free energy is additive, and thus the total free energy is the sum of the free energies of the 
subsystems. This line of reasoning is then implemented by the use of the replica technique, 
in an exact way for mean field models or in a variational way for more realistic systems 
The replica technique has been mostly tested for mean field models and references 
therein. 

Recently, a disorder dependent variational approach has been proposed for a large class of 
quenched disordered systems [f|. Roughly speaking, this method does the converse of replica 
variational methods: one looks for the best translationnally invariant Hamiltonian, for a 
fixed disorder configuration. One major advantage of this approach is that it is genuinely 
variational, thus providing a true upper bound to the free energy of the system, unlike 
replica based methods, which are "plagued by the n = limit". In fact, this method is 
closely related to the variational replica method, as shown for the random sine-Gordon || 
or the random-field XY model ||. In this paper, we apply this method to the problem of 
a d dimensional directed polymer chain in a random medium (or potential). This problem 
is related to many non-linear and/or disordered systems, such as random manifolds ||1D|| , 



kinetic growth ||11|| , turbulence ||12j| , and we hope that the present approach can shed some 



light on these problems. The Hamiltonian of the chain reads (in the following f3 — ^, where 
T is the temperature): 

m = 7^[ L ds 'r{sf + P f L ds V(f(s),s) (1) 
2a z Jo Jo 

where a is some microscopic length, f(s) denotes the d dimensional position of link s of the 
chain (s G [0,L]), and V(f(s),s) denotes the random potential acting on this link. In this 
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paper, we will restrict ourselves to a (site uncorrelated) Gaussian potential, the correlation 
function of which reads: 



V(r, s) V(P, s') = V 2 8(f- 7)5(s - s') (2) 

where A denotes the average of a quantity A over the disorder. 

The layout of this paper is the following. In section II, we define the variational Hamil- 
tonian, and calculate the corresponding variational free energy. Two different classes of 
variational solutions are exhibited. In section III, we discuss in detail the case of space 
dimension d < 2, where these two classes yield two different polymer width exponent v. In 
section IV, we study the case d > 2, , where we get an exponent v = | for all solutions. We 
argue that for d > 4, one of the variational solutions points undergoes a phase transition 
which is not unlike the (simple or generalized) random energy model. No evidence for a 
phase transition is found for 2 < d < 4 : one of the variational solutions can be interpreted 
in terms of Imry-Ma segments. In the conclusion, we discuss the implications of these re- 
sults, and suggest that a renormalization procedure should be performed together with this 
variational method. 



II. THE VARIATIONAL FREE ENERGY 

We consider Hamiltonian ([I]) and its associated variational Gaussian companion: 

(3H = lf L ds [ L ds' (f(s) - R) g-\s - s') (f(s') - R) (3) 
2 Jo Jo 

where we have restricted the variational kernel g to be translationally invariant. Note that 
contrary to the variational replica method where the disorder is a priori averaged, one has 
to consider here the possibility of a variational shift R for each link. (It would be even 
better to consider a variational link-dependent parameter R(s), but the calculation are too 
involved). The true free energy F satisfies the bound: 

f3F < p$(V) =(3F + (3<n-n >o (4) 
3 



where < ... >o stands for the thermal average with Hamiltonian (|j). Using equations ([!]), ([3|), 
and (Q), together with periodic boundary conditions r(0) = r(L), we obtain the disorder 
dependent variational free energy as: 

P*(V) = -dJ2 log 9 - 2 + d £ W ^ + /3W(-R, G) (5) 



n=l 1 ra=l 



with u> n = ^ and 



G) = / ^ f L ds V{k, s) e iM e-7 ^ (6) 



(27T) 







and 

00 00 / £ \ 

G = 2Efo = 2X; / da 0(5) e**" (7) 

n=l n=l V- 70 / 

Since .R is independent of s, and using equation (Q), we get 

W(R, G) = L h J-0T d U(k) e* 5 e-f ^ (8) 



with U{k)U{k') = V 2 5(k + k') 



A. The variational solutions 

The minimization equations with respect to g n and .R read: 

a = a 2 (9) 

W + ^V%W(R, G) 

and 

V n W(R,G) = (10) 

Note that one may find several solutions to the variational equations. In disordered systems, 
one does not expect, in general, that their (variational) free energy differ by an extensive 
amount. Since fluctuations around one solution, or instantons connecting different solutions, 
may yield extensive contributions to the free energy, one has a priori to keep all variational 
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solutions, unless some can be shown to be unstable with respect to such fluctuations 0. In 
the present problem, we first point out that there are two classes of solutions as is clear from 
equation (|9|): 

(i) the first, hereafter denoted by (+) has V|W( J R+, G+) > 0. 

(ii) the other, denoted by (-), has V|W(.R_,G_) < 0. 



Defining a + = |^-V|W(-R + , G + ), equation (|7|) can be rewritten 



27T 2 ^[ n 2 + a + 

for the (+) solution. In a similar way, we define a_ = 
(0) can be rewritten 



(11) 



f3La 2 

4ir 2 d I v R 



La 



2 oo 



E 



2tt 2 — [ n* - a_ 



V 2 ? W(^_, G-)\, and equation 



(12) 



for the (— ) solution. Note that we must have a_ < 1 for stability reasons, since equation 
(||) requires that g n > 0, Vn. Using equations (1.421.3 and 1.421.4) of reference jEJ, we 



may rewrite equations (O) and (1I2|) as 



(7_ = — ^ = (coth(7r v /ol 



47T -v/cEf 



7r v 'a + ' 



(13) 



and 



G_ 



La 2 1 1 

(cot [7T voT) 



47r 1 /oT 



The variational free energies 3>±(V) are easily obtained through equation 
#±(V) = $±(V) - $(0), we get 



(14) 



). Denoting 



{3m + = /3(W(R+, G+) - -G + V%W(R+, G+)) + log 



,sinh(7r v /tt 



) 



(15) 



for the (+) solution and 



p*_ = 0(W(iL, G_) + -G_|V|W( J R_, GL)|) + log 



sin 7r,/a_ 



(16) 



for the (—) solution. 



B. On the number of solutions and typical values of random quantities 



In principle, one should solve the minimization equations to find the disorder dependent 
quantities of interest G± and R±. In practice, we will estimate typical orders of magnitude 
in the following way. For instance equation M) implies 



/n ki r n a k -t -» n - 

(2^p J (2^p Uih) U{h) eikvRl+ik2 ' R2 e ~ f {kl+kl) ( 17 ) 

Clearly, we are only able to estimate this quantity if we temporarily forget that G and R 
depend themselves on U(k). This decoupling procedure then yields 
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LV (.R1-R2) 



W(R 1 ,G)W(R 2 , G) ~ ———i e — (18) 

{2^nG)2 

so that a typical value of W(R, G) reads 

{W(R,G)) t ^ V * (19) 
up to a random (algebraic) constant. In the same conditions, we obtain 

^v^n (20) 



Hyp ^ 



G A ' 4 



up to a random (positive) constant. 

To get a feeling for its range of validity, one may calculate, within this approximation, the 
averaged number of points R that satisfies the minimization equation (|10D, or equivalently 
the averaged (or typical) distance between such two such points R\ and R 2 . A straightfor- 
ward calculation, based on equations fllCf ), and (|18D shows that (\Ri — i? 2 |J ~ VG. 

We thus expect that the disorder dependent variational method, and the approximate 
estimation of typical random quantities are justified if there are few solutions, that is if G± 
is large. 



C. Stability of the solutions 

We will also examine the stability of the variational solutions with respect to the varia- 
tional parameters g n and R. As explained in reference for the random field XY model, we 
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do not expect any instability with respect to the g n 's and consider only the stability of the 
solution (G, R) with respect to small link-dependent shifts <5(s). The associated free energy 
reads 

A(/M>) =Jds (t^/(s) - /*/ (0? (* ' *Y V& s) e** e-f * j (21) 

The positivity of the A(/3$) is determined by the spectrum of the kernel in (|2"T|). This kernel 
is analogous to that of a Schrodinger equation in a random potential of typical strength 
( v d ) . For large values of G, the potential vanishes, leading to marginal (zero energy) 
fluctuations. 



III. RESULTS FOR D < 2 



For physical purposes, this means essentially d — 1. In this case, one has an exact 
solution |L4j] with v = | (and corrections to the free energy of order L^), and a variational 
replica calculation with v = | (and corrections to the free energy of order L^). We now 
consider the two classes of variational solutions, which (in d = 1) can be called potential 
minima or maxima. 



A. The (+) solution 



It is easily checked, using ([13]), that for d — 1, the only self consistent solution of equation 
(pT|) is a + ~ L§ and 

G + ~ L§ (22) 
that is an exponent = |, or more generally v = nzdj for d < 2. As shown in section 



( [11 CD , this solution is marginally stable. Moreover, its physical meaning can be appreciated 



through a Flory argument; since the (+) solution corresponds to attractive ("collapsed") 
regions, a balance between the entropic term and an Imry-Ma estimate [TB[ of the potential 
term yields 



which indeed yields v = j^zjj • For d = 1 , the disorder dependent part of the variational free 
energy \E' + is of order . 

B. The (-) solution 

From equation (0), one must have a_ < 1. This implies for d = 1 

G_ ~ Ll (24) 



yielding u = |, or more generally z/ = for d < 2. This solution is also marginally stable. 
It can also be obtained, in a Flory like manner, applied now to the repulsive (or swollen) 
regions 

In this case, the disorder dependent part of the variational free energy (see equation (JT6|) ) 
is of order Ls, in agreement with the Flory estimate (and the variational replica result). 

These results strongly suggest that the (— ) solution is very similar to the full symmetry 
breaking replica solution. The second length scale which comes out of the (+) (or domain) 
solution has not been obtained by the other methods. Note that its free energy is of 
order . Since G diverges with L for both solutions, we expect our variational approach 
to be meaningful: (i) there are few such solutions (ii) these solutions are marginally stable. 

IV. RESULTS FOR D > 2 

The problems we face for d > 2 are threefold: 

(a) in our approximations, the disorder becomes almost irrelevant for L large and G large. 



In marked contrast with the d < 2 case, equations ( |Il3| , |T4|) give a single solution, G ~ L, 
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(2-d) 

together with a± ~ VL 4 ~ 0. The fact that the exponent v sticks to its Brownian value 
above two dimensions has been also obtained in the variational replica method ||. 

(b) the identity of solutions (+) and (— ) does not survive if we allow for variational 
solutions where either G or L, or both, become finite. In this case, the variational method 
we have used requires at least a new interpretation, since it has many variational solutions 
(if G is finite), or considers only a finite portion of the chain (if L is finite), or both. It is 



also possible that such solutions are unstable (see section ( [II CD ). 

(c) to make matters worse, most of the high dimension models deal with directed poly- 
mers on a lattice. Most prominent among these lattice models are the three-dimensional 
and infinite-dimensional (tree) models ( [0] and references therein). It is clear that the com- 
parison of the continuous model described by equation ([[J) with these discrete models is not 
obvious, notwithstanding the very existence of the continuum limit jl 1] . 



With all these caveats in mind, we will now discuss two particular solutions of the 
variational equations, which may bridge the gap between discrete and continuum models. 



A. The large L, small G, (+) solution for d > 4 

It is possible to find a solution G + ~ for large L, as seen from equation (|l^) in the 

g 2 

limit of large a + . We get — £ ~ L< 4 ~ d ) which indeed vanishes for a > 4. (This solution 
corresponds, for d < 2, to the domain solution). Even though we do not wish to discuss 
in detail how the limits L large and a small are to be taken, we will see below that indeed 
d = 4 is a borderline dimension. 

The disorder dependent part of the variational free energy ( |IIf) reads for large a + 

~ W(R +1 G + ) - l -G + V\W{R + , G+) + i(V5T) (26) 

The first term on the r.h.s. of equation ([26|) is to be evaluated with the constraint that 
V^W(i? + , G + ) > 0. It is then easily shown that the algebraic constant omitted in equation 
(IBJ) is negative. The second term on the r.h.s. of (B^) is also negative. Finally, the third 



9 



term is positive: since we have = Q + (V) — $(0), this implies the existence of a phase 
transition between a high temperature Brownian phase — £ ~ L and a low temperature 

C 2 

"frozen phase" ~ L 4 - d . Note that the free energy of this frozen phase scale like y/di+, 
i.e. like L< d - 4 ) , which explicitly shows the problems of the continuum version of the model 
as compared to its lattice counterpart. 

As stressed above, this phase transition pertains to a single (+) solution. The typical 
distance between two "frozen" solutions being of the order of a, we are thus faced with an 
exponential number (in L) of such solutions. To get a flavour of the nature of the phase 
transition, we calculate a typical correlation between the free energies of two such solutions, 
and obtain, within the decoupling scheme of section ( [II B| ): 

- - LV 2 ( ^i-R2) 2 \ , 

W(R 1 ,G + )W(R 2 ,G + )^ * el 4G + ) ~ n^LV 2 5(Ri — R 2 ) (27) 

(2 x /7rG + )2 

These correlations are indeed reminiscent of the (simple or generalized) random energy 
model. We cannot evaluate the exponent v in the low temperature phase, but we think it is 
also v — ^, since the polymer undergoes a random walk between variational frozen solutions, 
that is a random walk on the different (R4) points. We tentatively conclude that dimension 
d = 4 may well be some kind of lower critical dimension for a (simple or generalized) random 
energy model type of phase transition. The role of dimension d = 4 in this context has been 
recently discussed in |i~6| . 



B. The finite L, finite G, (— ) solution for d > 2 

Another intriguing result of the variational equations concerns the (— ) solution, since it 
is restricted by the condition a_ < 1. This condition does not play any role in the (large 
L, large GL) solution, for d > 2. However, if one considers, a finite G_ solution, it can only 



exist up to a number of links L < Lq, (see equation (I2[) and the definition of a_) with 
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4vr 2 d ( f d d k ^ t ^ fe 2 N ~' 



P \J (2tt) 



= — [ I 7^-T d k2 U(k) e lk - R - e-T - ) (28) 
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Knowing the probability distribution of U(k), and using the decoupling approximation of 
section ( [11 B| ), one may evaluate the probability distribution of L as 



1 4- — \ / 24- — 

G \ I G 



P(Io) .Mil. IV) j-^-j exp-^C(d,/3y)-^-j (20) 

where A(d, (3V) and C(d, f3V) are regular functions of the dimension and of the temperature. 
Note that the second (and higher) moment of this distribution is divergent. This result 
strongly suggests that the chain breaks into Imry-Ma domains of (distributed) size Lq. A 
more detailed description (such as the role of the scale L in the overlap between different 
variational solutions) requires a better understanding of the spatial succession of (+) and 
(— ) finite L solutions along the chain. 

V. CONCLUSION 

In this paper, we have presented a disorder dependent variational method for the problem 
of a d dimensional directed polymer in a random potential. This method seems reliable for 
d < 2 , where there are few variational solutions, and agrees when they overlap, with the 
(full) replica symmetry breaking variational method. We have also found a new length scale 
(the domain solution), which is apparently missed by other approaches. For d > 2, the 
variational solutions may be very dense and our variational procedure should be viewed as 
a first step towards a variational renormalization group: the free energy ^(R, G) has indeed 
(see equations (|l^), ([16|) ) the form of a new random potential, so one may think of iterating 



the process ||17|| . In this approach, we have presented some peculiar solutions which may 
have some relevance, either to the puzzle of critical dimensions for this problem, or to the 
physical description of the chain. We have explicitly shown that d = 4 plays a special role 
for the (+) solutions, and that an Imry-Ma length L Q is, for d > 2, the natural scale for 
the correlations between different finite G_ solutions. A more ambitious goal would be to 
study the spatial interplay of the (+) and (— ) solutions in the variational renormalisation 
procedure to see if chaotic behaviour may arise [IB]. 
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